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Abstract  Understanding  the  adsorption  of  water  vapor  onto 
activated  carbons  is  important  for  designing  processes  to  re¬ 
move  dilute  contaminants  from  humid  gas  streams,  such 
as  providing  protection  against  chemical  warfare  agents 
(CWAs),  or  against  toxic  industrial  compounds  (TICs) 
used  in  a  terrorist  chemical  attack.  Water  vapor  isotherms 
for  Calgon  BPL  granular  activated  carbon  (GAC),  mil¬ 
itary  ASZM-TEDA  GAC,  electrospun  activated  carbon 
nanofibers  (ACnF),  Calgon  Zorflex™  activated  carbon  cloth, 
and  Novoloid-based  activated  carbon  fiber  cloth  (ACFC)  are 
presented.  Of  particular  interest  are  the  ACFC  isotherms, 
which  exhibit  an  unusually  high  degree  of  hydrophobicity. 
The  ACFC  isotherms  also  show  a  correlation  between  water 
vapor  adsorption  hysteresis  and  the  level  of  activation.  Water 
vapor  isotherm  models  from  the  literature  are  compared. 

Keywords  Activated  Carbon  Fiber  Cloth  (ACFC)  •  Water 
vapor  adsorption  •  Hysteresis  •  Isotherm  models 

Abbreviations 

A  Distribution  constant  for  Mahle  equation 

a  Amount  adsorbed  at  a  given  pr 
ac  Constant  for  DS4  equation 

a  f  Capacity  of  functional  site  part  of  Do  &  Do  type 
equations 

a:  Capacity  of  type  I  term  of  Stoeckli  formula,  or 

capacity  of  Ising  term  in  Extended  CMMS  equation 
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ai  Capacity  of  Langmuir  term  of  Extended  CMMS 
formula 

am  Capacity  constant  in  Talu-Meunier  equation 
a0  Pore  capacity 

ap  Concentration  of  primary  sites 

ay  Capacity  of  type  V  term  in  Stoeckli  equation 
Capacity  of  microporous  part  of  Do  &  Do  type 
equations 

B  Distribution  constant  for  Mahle  equation 
c  Henry’s  Law  coefficient  for  DS  equations 
E0  Reference  energy 

H  Proportionality  constant  in  Talu-Meunier  equation 

K  Constant  for  Lagorsse  equation 

Ka  Langumuir-type  constant  for  adsorption  part  of 

Lagorsse  equation 

Kas  Shape  constant  for  CMMS  equation 
Kd  Langmuir-type  constant  for  desorption  part  of 
Lagorsse  equation 

K  f  Langmuir-type  constant  for  functional  site  part  of 
Do  &  Do  type  equations 

K/jL  Langmuir  type  constant  for  microporous  part  of  Do 
&  Do  type  equations 
k  Constant  for  QHR  equation 

kQ  Langmuir-type  constant  for  CMMS  equations 

k\  Constant  for  CMMS  equations 

kL  Langmuir  constant 

kt  Constant  in  Talu-Meunier  equation 

N  Number  of  terms  to  sum 

n  Dubinin-Astakhov  exponent 

ni,ny  Stoeckli  exponents 
p$o  Relative  pressure  at  50%  of  a0  adsorbed 
pr  Relative  pressure  of  water  (relative  humidity) 

Rg  Universal  gas  constant 

T  Temperature 
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P  Energy  coefficient  of  an  adsorbate-adsorbent 

system  (fiE0  =  E,  the  affinity  energy  of  the  system) 
Pi,  Pv  Energy  coefficient  of  each  term  in  the  Stoeckli 
formula 

k  Fitting  parameter  for  Dubinin-Serpinsky  type 
equations 

Fitting  parameters  for  LeVan  equation 

1  Introduction 

Understanding  water  vapor  adsorption  onto  activated  car¬ 
bons  is  important  for  adsorption  processes  such  as  steam 
regeneration  of  activated  carbons  or  the  removal  of  contami¬ 
nants  from  humid  gas  streams.  Water  can  compete  for  avail¬ 
able  adsorption  sites  and  reduce  the  capacity  for  the  carbon 
to  capture  target  compounds.  Under  high  humidity,  water 
can  completely  fill  the  micropores,  nearly  eliminating  the 
capacity  to  capture  any  contaminants  more  weakly  adsorbed 
than  water  (Huggahalli  and  Fair  1996).  The  coadsorption  of 
water  and  organics  onto  GAC  has  been  studied,  including 
the  hysteresis  effect  (Rudisill  et  al.  1992). 

One  relevant  application  for  which  contaminants  are  re¬ 
moved  by  adsorption  from  humid  gas  streams  is  providing 
protection  against  a  chemical  attack  in  either  military  or 
homeland  defense  (terrorist)  situations.  Whether  purifying 
air  with  a  respirator  for  personal  protection,  or  purifying  air 
with  a  filtration  system  for  collective  protection  of  a  shelter 
or  building,  outside  air  is  the  medium  to  be  filtered,  and  it 
is  not  humidity-controlled.  Another  application  where  out¬ 
side  air  may  be  filtered  is  in  large-scale  painting  operations, 
where  exhaust  volumes  can  range  to  the  thousands  of  cu¬ 
bic  meters  per  second,  and  the  cost  of  humidity  control  is 
prohibitive. 

Activated  carbon  fibers  (ACFs)  show  great  promise  as 
a  new  adsorbent  and  catalyst  support  in  competition  with 
more-traditional  activated  carbons  (Mays  1999).  Carbon 
fibers  can  be  woven  into  cloth  and  then  activated  to  form 
activated-carbon-fiber  cloth  (ACFC).  ACFCs  are  generally 
made  from  pitch,  rayon,  polyacrylonitrile  (PAN),  or  pheno¬ 
lic  (Novoloid,  Hayes  1981)  resin.  The  fibers  are  woven  into 
a  fabric  and  then  activated  in  steam  or  CO2  at  high  temper¬ 
atures  (e.g.,  800  °C).  ACFCs  have  the  ability  to  be  regener¬ 
ated  by  direct  electrical  resistance  heating.  Resistance  heat¬ 
ing  has  been  shown  to  be  a  more  rapid  and  energy-efficient 
method  of  regenerating  activated  carbon  than  desorbing  with 
steam  or  inert  gas  (Petkovska  and  Mitrovic  1994).  Because 
of  this,  ACFCs  have  been  studied  recently  in  an  effort  to 
develop  next-generation  energy-efficient  cyclic  adsorption 
systems  (Subrenat  and  Fe  Cloirec  2006;  Vidal  et  al.  2006; 
Hashisho  et  al.  2007). 

The  reduced  organic  adsorption  in  humid  gas  streams 
onto  ACFCs  has  been  studied  for  both  equilibrium  (Cal  et 


al.  1996)  and  dynamic  (Sullivan  et  al.  2001)  adsorption.  Wa¬ 
ter  adsorption  isotherms  are  herein  presented  for  a  variety 
of  ACFCs,  and  Calgon  BPF  as  a  reference,  in  an  attempt 
to  more  thoroughly  understand  the  adsorption  of  water  onto 
these  materials. 

2  Experimental  methodology 

2.1  Carbon  selection  and  sample  preparation 

Carbons  examined  in  this  study  included  Calgon  BPF™ 
granular  activated  carbon  (GAC),  military  ASZM-TEDA 
granular  activated  carbon  (GAC)  from  a  C2A1  respirator, 
Activated  carbon  nanofiber  (ACnF),  Calgon  Zorflex™  ac¬ 
tivated  carbon  cloth,  and  Novoloid-based  activated  carbon 
fiber  cloth  (ACFC).  The  matrix  of  Novoloid-based  ACFCs 
tested  consisted  of  four  levels  of  activation  of  virgin,  oxi¬ 
dized,  and  H2 -treated  fabrics. 

The  ACFC  used  was  American  Kynol  ACC-5092- 10, 
ACC-5092-15,  ACC-5092-20,  and  ACC-5092-25.  The  last 
two  digits  of  the  ACFC  part  number  are  related  to  the  level 
of  activation  and  are  analogous  to  the  nominal  surface  area 
of  the  material,  i.e.,  -10  is  ~1000  m2/g,  -15  is  ~1500  m2/g, 
and  so  on.  The  Calgon  Zorflex™  material  tested  was  ACC- 
FM10  and  ACC-FM 10-700.  The  ACnF  material  used  was 
produced  by  eSpin  Technologies  Inc.,  of  Chattanooga,  Tenn. 
The  ACnF  is  a  non-woven,  produced  by  electrospinning 
polyacrylonitrile  (PAN),  and  then  activated  by  conventional 
means. 

The  oxidized  samples  were  obtained  by  treating  virgin 
ACFC  samples  (0.025  m2,  4  to  7  g)  having  different  grades 
of  activation  with  150  mL  of  1/1  (V/V)  concentrated  HNO3 
and  H2SO4  solution  for  5  days  (Dimotakis  et  al.  1995a, 
1995b).  Samples  were  then  washed  with  distilled  water  un¬ 
til  neutral  pH  conditions.  Samples  were  then  drained  from 
water  and  dried  in  an  oven  at  150  °C  in  a  flow  of  200  seem 
of  N2. 

The  H2-treated  ACFC  samples  were  heated  at  950  °C 
in  pure  hydrogen  environment  (Menendez  et  al.  1996; 
Fi  et  al.  2002).  Each  sample  was  put  in  a  quartz  boat  inside  a 
5 -cm  outer-diameter  quartz  tube  fitted  inside  a  temperature- 
controlled  tube  furnace  (Findberg  furnace  model  54232; 
temperature  controller  unit  59344).  The  samples  were 
heated  in  N2  for  2  hours  at  200  °C  to  remove  chemisorbed 
oxygen  and  then  the  gas  was  switched  to  hydrogen  (Zero  gas 
grade)  at  200  seem.  The  samples  remained  at  950  °C  for  3 
hours,  then  heating  was  stopped  and  the  samples  were  left  to 
cool  down.  When  the  temperature  decreased  below  100  °C, 
the  blanketing  gas  was  switched  to  N2. 

Samples  were  labeled  according  to  the  precursor  fol¬ 
lowed  by  a  suffix  referring  to  the  type  of  treatment,  for 
instance  -15V  is  untreated  (virgin)  ACC-5092-15,  -15A  is 
acid-treated  (oxidized)  ACC-5092-15,  -15H  is  hydrogen- 
treated  ACC-5092- 15. 
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2.1.1  N2  isotherms,  surface  area  and  pore  size  distribution 

N2  adsorption  isotherms  were  generated  by  volumetric  mea¬ 
surements  using  ASAP  2010  instrument  (Micromeritics, 
Inc.).  Samples  were  degassed  for  more  than  8  hours  at 
150  °C  under  vacuum  until  a  pressure  of  4-6  pm  Hg.  Acid- 
treated  samples  have  to  be  degassed  for  longer  periods  com¬ 
pared  to  the  virgin  or  the  hydrogen-treated  samples  to  reach 
this  degassing  pressure  limit.  The  adsorption  isotherms  were 
obtained  at  a  relative  pressure  of  10-6  to  0.99.  N2  adsorption 
and  desorption  occurred  at  the  N2  boiling  point  (77  K). 

The  standard  Brunauer-Emmett-Teller  (BET)  method 
was  used  to  fit  the  adsorption  data  over  a  relative  pressure 
range  of  0.06  to  0.20.  Micropore  volume  and  pore  size  dis¬ 
tribution  were  determined  using  the  3D  model  (Sun  et  al. 
1998).  The  model  is  based  on  the  Dubinin-Radushkevich 
(DR)  equation  and  is  used  to  model  the  adsorption  of  gases 
on  microporous  adsorbents  with  slit  shaped  pores  between 
4  A  and  100  A.  The  average  micropore  width  was  deter¬ 
mined  based  on  pore  size  distribution  for  sizes  <20  A  while 
the  total  pore  volume  was  determined  based  on  the  volume 
of  nitrogen  adsorbed  at  saturation  and  using  nitrogen  liquid 
density  (0.808  cm3/g). 


2.1.2  X-ray  photoelectron  spectroscopy  (XPS) 

The  elemental  content  (C,  O,  N,  and  S)  at  the  surface  of 
the  virgin  and  chemically  treated  ACFC  samples  was  de¬ 
termined  using  XPS  technique.  The  analysis  was  conducted 
using  a  PHI  5400  (Perkin-Elmer,  Physical  Electronics  Inc.) 
instrument.  A  spot  size  of  1  mm2  was  targeted  with  Mg  Ka 
radiation  or  Al  Ka  radiations  at  a  power  of  300  W  at  15  kV 
under  ultrahigh  vacuum  (10-8  to  10-9  torr). 

2.1.3  H2O  adsorption  isotherms 

H2O  adsorption/desorption  isotherms  were  collected  gravi- 
metrically  with  a  Gravimetric  Sorption  Analyzer,  Model 
GHP-F,  manufactured  by  VTI  Corporation  of  Hialeah,  Fla. 
(Fig.  1).  The  system  uses  a  Magnetic  Suspension  Balance  by 
Rubotherm  of  Bochum,  Germany.  The  matched  mass  flow 
controllers  in  the  system  are  rated  for  a  nominal  1.0  sFpm 
flow.  Carbon  samples  used  for  the  isotherms  ranged  from  40 
mg  to  150  mg. 

Nitrogen  was  used  as  the  carrier  gas  throughout  the  ex¬ 
periments.  Each  sample  (except  the  ASZM-TEDA  GAC) 
was  dried  at  250  °C  under  nitrogen  prior  to  each  isotherm 
until  a  static  weight  was  achieved  (<0.03%  change  in  10 


Fig.  1  Gravimetric  sorption  analyzer  model  GPH-F  schematic  and  photo 
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minutes)  or  for  a  maximum  time  of  one  hour.  The  ASZM- 
TEDA  carbon  was  dried  at  105  °C  to  minimize  impact  on 
the  impregnants  in  the  carbon. 

The  GHP-F  contains  algorithms  for  generating  a  spec¬ 
ified  relative  humidity  (RH)  based  on  the  evaporator  tem¬ 
perature,  the  sample  temperature,  the  flow  rates  of  the  mass 
flow  controllers,  and  a  correlation  for  the  vapor  pressure  of 
water.  To  verify  the  accuracy  of  the  RH  generated  by  the  in¬ 
strument,  an  Edgetech  Vigilant™  Chilled-Mirror  Hygrome¬ 
ter,  with  a  model  S3  sensor,  was  connected  in  series  with  the 
GHP-F.  The  set-point  and  measured  RH  values  were  within 
an  average  1%  RH.  The  total  gas  flow  through  the  system 
was  monitored  with  a  mass  flow  meter,  Omega  FVL-1600 
0-5  LPM. 


3  Equilibrium  isotherm  modeling 

A  large  collection  of  mathematical  formulas  have  been  pro¬ 
posed  by  previous  investigators  to  model  isotherms  of  ad¬ 
sorption.  Here,  we  present  an  overview  of  equations  that 
have  been  proposed  to  model  adsorption  of  water  onto  car¬ 
bon.  Many  of  the  more  recently  developed  models  have  been 
fit  to  only  a  few  sets  of  experimental  data.  While  our  ap¬ 
proach  to  these  formulas  is  genealogical,  we  will  fit  all  these 
equations  to  our  data  to  examine  which  equations  better 
model  the  spectrum  of  materials. 

We  will  attempt  to  elucidate  the  evolution  from  one  for¬ 
mula  to  another  and  briefly  state  similarities  of  form  to  other 
equations.  To  make  these  similarities  more  clear,  we  have 
sometimes  stated  the  equations  in  a  different  form  from 
those  given  by  the  original  authors.  Quite  a  few  recent  equa¬ 
tions  take  the  form  of  a  summation  of  two  terms;  these 
include  the  Stoeckli  equation,  all  equations  in  the  Do  & 
Do  family,  and  the  extended  CMMS  equation.  Oftentimes, 
when  fit  to  a  set  of  experimental  data,  one  of  the  factors  pre¬ 
ceding  the  term  representing  how  much  of  the  adsorption 
that  term  is  responsible  for  (represented  by  a  with  a  sub¬ 
script)  will  become  zero  or  close  to  zero;  this  indicates  that 
the  term’s  contribution  is  insignificant,  and  overfitting  is  oc¬ 
curring.  Usually  if  this  factor  is  set  to  zero,  the  equation  will 
reduce  to  a  previously  developed  isotherm  formula. 

It  is  important  to  evaluate  these  models  not  only  on  a 
quantitative  basis — that  of  measures  of  error  and  correla¬ 
tion — but  also  a  qualitative  basis.  Mathematical  qualities 
that  should  be  satisfied  including  having  a  finite  (usu¬ 
ally  nonzero)  Henry’s  Law  limit,  some  certain  fixed  val¬ 
ues  (a  =  0  at  pr  =  0,  a  =  aQ  at  pr  =  1),  and  increasing 
monotonically  over  the  entire  range  of  RH  values.  They 
should  also  be  independent  of  units  used.  It  is  also  worth 
mentioning  that  some  formulas  may  fit  to  experimental  data 
better  more  due  to  the  flexibility  of  the  form  of  their  equation 
than  to  the  theory  used  to  derive  them.  If  a  formula  is  fit  to 


a  data  set  from  a  carbon  with  properties  very  different  from 
the  assumptions  made  in  deriving  the  model,  it  may  produce 
a  good  fit  regardless,  though  the  fitting  parameters  may  not 
have  physical  meaning.  This  is  increasingly  common  as  the 
number  of  parameters  to  a  formula  becomes  larger.  Since 
many  proposed  formulas  are  modifications  to  previous  for¬ 
mulas,  we  will  group  them  and  examine  them  as  such. 

In  Table  1,  pr  is  the  relative  pressure  of  water  (relative 
humidity),  a  is  the  amount  adsorbed  (in  mg/g),  T  is  tem¬ 
perature  (in  Kelvin),  Rg  is  the  universal  gas  constant,  E0  is 
a  reference  energy,  and  A  is  a  fixed  value  indicating  how 
many  terms  to  sum  over.  All  other  values,  unless  defined 
below,  are  treated  as  fitting  parameters;  except  for  £z-  in  the 
LeVan  equation,  all  fitting  parameters  must  be  positive. 

4  Dubinin-Radushkevich  type  formulas 

The  term  ‘micropore  filling’  was  first  coined  in  Dubinin’s 
adsorption  work;  since  micropore  filling  is  important  to 
these  carbons,  it  is  reasonable  to  start  there.  The  origi¬ 
nal  Dubinin-Radushkevich  equation  dates  to  1947.  The  DR 
equation  was  found  to  perform  well  on  most  activated  car¬ 
bon,  but  not  to  deal  well  with  carbons  with  very  fine  micro¬ 
pores  (Do  1998,  p.  156).  The  Dubinin-Astakhov  equation 
replaces  the  exponent  of  2  with  n ,  a  parameter  that  increases 
with  the  degree  of  activation. 

The  Stoeckli  equation  (Stoeckli  2002),  created  to  fit  type 
IV  isotherms  on  oxidized  carbon,  is  a  sum  of  two  DA  equa¬ 
tions,  one  to  account  for  the  type  I  behavior  and  one  to  ac¬ 
count  for  the  type  V  behavior. 

The  Doong-Yang  equation  (Doong  and  Yang  1987)  is 
one  of  only  two  equations  we  have  found  that  accounts  for 
hysteresis  in  desorption.  To  model  the  desorption,  hQ  is  set  to 
the  relative  pressure  at  the  beginning  of  the  hysteresis  loop. 
Gamma  is  an  empirical  activity  coefficient,  and  is  found  us¬ 
ing  the  Dubinin-Serpinsky  equation  DS2,  which  will  be  dis¬ 
cussed  in  the  next  section. 

All  equations  of  this  family  have  a  Henry’s  Law  limit  of  0 
(as  long  as  n  >  1)  (Do  1998,  p.  162),  which  many  previous 
investigators  have  remarked  is  thermodynamically  inconsis¬ 
tent.  However,  in  very  hydrophobic  carbons,  the  value  of  the 
Henry’s  Law  constant  is  very  small,  and  this  limit  can  be  a 
reasonable  approximation. 

5  Dubinin-Serpinsky  type  formulas 

The  first  four  equations  of  DS  type  vary  in  the  factor  that 
is  used  to  account  for  the  decrease  of  available  adsorption 
sites.  DS2  was  proposed  by  Dubinin  (Dubinin  and  Serpin- 
sky  1981),  and  DS4  by  Barton  et  al.  as  a  modification  to 
the  DS2  (Barton  et  al.  1992).  There  are  members  of  the  DS 
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Dubinin-Radushkevich  type  formulas 
Dubinin-Radushkevich  (DR) 
Dubinin-Astakhov  (DA) 

Stoeckli 

Doong-Yang  (DY) 
Dubinin-Serpinsky  type 
Dubinin-Serpinsky  (DS2) 

DS4 

LeVan 

Do  &  Do  type 
Do  &  Do 
(DD) 

CIMF 

Cossarutto 

Lagorsse 


Other 

Qi-Hay-Rood 

(QHR) 

Mahle 

Talu-Meunier 

CMMS 


CMMS  Type  V 
(Ising) 

Extended  CMMS 
(Ising-Langmuir) 


a  =  acexp[-(^)2],  A  =  RgT ln^ 
a  =a0  exp[—  ( )"  ] 

a  =  a7  exp[-(^)n']  +  av  exp[-(^^)”v] 
a  =  acexp(-(|^)2[(ln  ^)2  -  (In ha)2]) 

Pr  ~  c(ap+a)(l-Ka) 

p  —  _ “ _ 

yr  cap+ca(l  — exp(— /c2(a— ac)2)) 


l+Kr  Y,i=i  pi 


K f  Yf=l  iPr 
a~af  l+KfYT+iPi 


11  \+Ktljf 
KnP? 


a-  nr  Kfpr  I  a  Wr 

a  — Of  i+KfPr  -T  “p,  1 +Kppnr 


a 


1  ELi M 


Ka  if 


Adsorbing:  a  =  afT^^^-+a^ 


Desorbing:  a  =  ay 


i  Ek  -W)  i  n  KdPJ 
1  +KdpJ  l+T.Ui(Kp>)  M  1 +KdpJ 


l+exp[k(p5o~Pr)] 


arctan(  PrB  A )— arctan(  -ff ) 
°  arctan  ( ) — arctan  ( ) 

i*/a„ 

ko£r 


Pr  =  J W  V  =  -1+™,  ?  =  ^ 


(l-Kas)(koPr+W2(l-Kas)) 

ft 


W  —  1(1  nf 


k\Pr 

(1  ~Kas) 


+  i 


k]_Pr 

(1  ~Kas) 


)2  + 


4k0pr 
(1  ~Kas) 


a=ao  .  U)=  5(1  -felPr+x/(T  -fclPr)2 +4*0pr) 


a=«LT^;+«/ 


2 

koPr 

koPr+W 2 


,  w  as  used  in  the  Ising  equation 


2nd  order 
3rd  order 
6th  order 

4th  order 

3rd  order 
4th  order 
Nth  order 

4th  order 

5  th  order 
5  th  order 
5  th  order 


3rd  order 

3rd  order 
3rd  order 

4th  order 
3rd  order 

5th  order 


family  we  have  chosen  not  to  treat.  DS1  (Dubinin  1980) 
fits  poorly  to  data  at  middle-to-high  RH  and  is  not  treated 
here.  DS3  (Barton  et  al.  1991)  claims  that  the  relative  hu¬ 
midity  when  a  =  0  is  usually  more  than  1/3,  and  that  the 
derivative  da /dp  there  is  large;  this  is  clearly  aphysical.  The 
Corrected  DS  equation  (CDS)  (Gauden  2005)  was  proposed 
by  Gauden  and  contains  an  extended  empirical  factor  to  im¬ 
prove  on  DS4  at  relative  pressures  less  than  10%,  but  due  to 
the  limitations  of  our  experimental  system,  we  do  not  have 
adsorption  data  at  relative  pressures  between  0%  and  10%. 
The  CDS  can  also  be  shown  to  not  be  independent  of  units, 
because  of  logarithms  and  exponentials  in  its  formula.  For 
these  reasons  we  have  chosen  not  to  treat  these  equations. 

The  LeVan  equation  (Qi  and  LeVan  2005)  is  a  very  flex¬ 
ible  formula  that  is  essentially  a  truncated  power  series  of 
n/p  in  terms  of  n.  It  is  included  here  because  the  derivation 


follows  that  of  the  DS  equations,  and  it  contains  the  DS2 
and  DS3  equations  as  special  cases  (N  =  2  and  N  =  3,  re¬ 
spectively).  It  is  also  evident  that  if  N  was  allowed  to  tend 
to  infinity,  it  would  be  able  to  contain  the  DS4  and  the  CDS 
equations,  as  well  as,  in  theory,  any  other  isotherm  formula. 
Usually  N  =  4  is  enough  to  get  a  good  fit.  Aside  from 
which  is  the  Henry’s  Law  constant,  the  physical  meaning  of 
other  parameters  to  the  equation  is  not  clear.  Also,  because 
of  its  flexibility,  care  must  be  taken  to  make  sure  that  its  pa¬ 
rameters  are  such  that  the  equation  is  increasing  and  positive 
throughout  its  domain. 

The  stated  form  of  these  equations,  which  give  the  rela¬ 
tive  humidity  in  terms  of  the  loading,  is  the  opposite  of  the 
usually  desired  form  and  is  awkward  to  use  in  practice.  DS2 
can  be  stated  to  give  the  loading  in  terms  of  the  RH  by  solv¬ 
ing  a  quadratic,  but  this  form  is  somewhat  unwieldy;  since 
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for  purposes  of  fitting  the  two  equations  are  equally  useful, 
we  have  chosen  to  use  the  stated  form. 


6  Do  &  Do  type  formulas 

The  equation  of  Do  and  Do  (2000)  to  model  water  isotherms 
has  received  a  good  deal  of  attention  since  its  proposal.  One 
of  Do  &  Do’s  key  assumptions  is  that  water  inside  micro¬ 
pores  forms  clusters  of  molecules.  The  original  Do  &  Do 
(DD)  formula  assumes  clusters  of  5.  The  CIMF  isotherm, 
proposed  shortly  afterwards,  generalized  Do’s  equation  to 
clusters  of  arbitrary  size,  and  Cossarutto’s  modification  to 
the  DD  formula  went  one  step  further  and  simplified  the  por¬ 
tion  of  the  equation  representing  adsorption  on  functional 
groups.  The  first  term  in  these  three  models  is  a  SIPS  equa¬ 
tion,  and  the  second  term  of  Cossarutto’s  equation  is  a  Lang¬ 
muir  isotherm. 

Lagorsse  et  al.  proposed  the  second  formula  that  accounts 
for  desorption  (Lagorsse  et  al.  2005).  They  follow  a  simi¬ 
lar  derivation  to  the  DD  equation,  but  assume  clusters  of  7 
molecules.  The  first  term  of  both  equations  of  the  Lagorsse 
formula  is  also  a  SIPS  isotherm. 

Furmaniak  et  al.  proposed  a  “Heterogeneous  Do-Do 
Model”  that  accounts  for  different  energies  of  adsorption  at 
different  functional  sites  (Furmaniak  et  al.  2005).  However, 
without  prior  knowledge  of  the  nature  of  functional  sites  on 
the  material,  one  must  make  guesses  as  to  the  number  and 
properties  of  the  sites  the  material  may  have.  Since  we,  un¬ 
fortunately,  do  not  have  that  knowledge,  we  have  chosen  to 
not  treat  this  model  in  our  paper. 

We  note  that  the  summands  in  these  equations  have  ex¬ 
plicit  forms  that  simplify  computation. 

7  Other  formulas 

The  Qi-Hay-Rood  formula  is  another  type  5  model,  derived 
to  be  an  improvement  on  the  DS4  equation.  The  QHR  is 
nonzero  at  pr  =  0,  but  its  value  there  is  very  close  to  zero 
(Qi  et  al.  1998). 

Mahle  proposes  a  type  5  formula,  based  on  the  assump¬ 
tion  that  the  pore  distribution  approximates  a  Cauchy  dis¬ 
tribution,  which  is  similar  in  shape  to  a  normal  distribution 
(Mahle  2002).  It  is  known  that  the  pore  distribution  in  acti¬ 
vated  carbons  is  often  bimodal,  so  this  assumption  may  not 
always  hold.  The  general  method  for  deriving  an  isotherm 
from  a  distribution  of  pore  sizes  is  discussed  in  chapter  3.10 
of  Adsorption  Analysis:  Equilibria  and  Kinetics  (Do  1998). 

Talu  and  Meunier  derived  a  type  V  isotherm  that  assumes 
water  molecules  are  more  attracted  to  each  other  than  to  the 
carbon,  and  therefore  once  water  is  adsorbed  to  a  functional 
group  other  waters  will  cluster  onto  that  original  water  mole¬ 
cule  (Talu  1996). 
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Fig.  2  Water  vapor  isotherms  with  the  SPC  model  for  slit  pores  of 
0.65  nm  (Liu  and  Monson  2006).  Reproduced  with  permission 

Rutherford  proposed  a  model  based  on  cooperative  multi- 
molecular  sorption  theory  (CMMS)  that  explains  well  the 
effects  of  interaction  with  functional  sites,  but  does  not  ac¬ 
count  for  microporosity  or  condensation,  and  does  not  work 
well  at  relative  humidities  above  80%  (Rutherford  2003). 
The  CMMS  model  can  account  for  type  II  and  type  V 
isotherms.  Assuming  a  type  V  isotherm,  the  CMMS  model 
reduces  to  an  equation  called  the  Ising  equation.  Rutherford 
also  proposed  an  Extended  CMMS  model  to  better  examine 
carbons  with  minimal  micro-  or  nanoporosity  (Rutherford 
2006).  The  form  of  the  extended  CMMS  equation  is  a  sum 
of  an  Ising  equation  and  a  Langmuir  equation. 

8  Molecular  simulations  of  adsorption 

With  increasing  amounts  of  computational  power  available 
to  the  researcher,  molecular  simulations  of  adsorption  have 
become  more  popular.  Earlier  work  is  summarized  by  Bren¬ 
nan  (Brennan  et  al.  2001).  More  recently,  Liu  &  Monson  us¬ 
ing  the  Gibbs  Ensemble  Monte  Carlo  (GEMC)  method  pro¬ 
duced  isotherms  for  a  carbon  assumed  to  consist  of  slit  pores 
with  no  active  sites,  and  postulated  that  a  microporous,  hy¬ 
drophobic  carbon  would  resemble  their  isotherms  (Liu  and 
Monson  2005).  We  will  later  show  that  the  isotherms  for  our 
hydrogen-treated  Kynol™  material  bear  certain  similarities 
to  these  idealized  molecular  simulations. 

In  a  more  recent  paper,  Liu  and  Monson  performed 
Monte  Carlo  simulations  on  a  carbon  made  up  of  graphite 
platelets  with  randomly  distributed  active  sites.  The  iso¬ 
therms  from  those  simulations  strongly  resemble  the  iso- 
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therms  of  Calgon  BPL  and  similar  carbons  (Liu  and  Monson 
2006). 


9  Results  and  discussion 

9.1  XPS  of  virgin  and  chemically  treated  ACFC 

Oxidizing  the  ACF  enriches  the  ACFC  with  oxygen-contain¬ 
ing  functional  groups  such  as  carboxylic  acids,  quinones, 
and  phenolic  hydroxides.  Treating  the  ACFC  samples  with 
the  mixture  of  concentrated  HNO3  and  H2SO4  increased 
the  average  oxygen  content  among  the  samples  from  4.5% 
to  23.9%  while  the  average  carbon  content  decreased  from 
94.7%  to  74.6%.  The  average  nitrogen  content  increased 
from  0.5%  to  1.2%  while  the  average  sulfur  content  in¬ 
creased  from  0.21  to  0.25%.  Since  the  change  in  the  sul¬ 
fur  and  nitrogen  content  of  the  oxidized  samples  was  small 
and  disproportionate  to  the  increase  in  the  oxygen  content, 
this  indicates  that  the  increase  in  oxygen  content  is  not  due 
to  sulfate  or  nitrate  groups  from  H2SO4  or  HNO3.  On  the 
other  hand,  heating  the  ACFC  samples  at  950  °C  in  UHP 
H2  decreased  the  average  oxygen  content  among  the  sam¬ 
ples  from  4.5%  to  1.5%,  the  average  nitrogen  content  from 
0.5%  to  0.4%  while  the  average  carbon  content  increased 
from  94.7%  to  98.1%.  For  acid  and  hydrogen  treatments, 
the  change  in  oxygen  content  didn’t  seem  to  depend  on  the 
activation  level  of  the  virgin  ACFC  samples. 

9.2  N2  adsorption,  surface  area,  and  pore  size  distribution 
of  virgin  and  chemically  treated  ACFC 

The  H-treated  ACFC  depicted  comparable  BET  surface  area 
to  virgin  ACFC.  Going  from  ACFC  10  to  ACFC25,  the  BET 
surface  area  increased  for  both  the  virgin  and  the  H-treated 
ACFC.  On  the  other  hand,  BET  area  was  comparable  for  all 
the  oxidized  ACFC,  and  was  lower  than  that  of  virgin  ACFC. 
The  BET  area  for  the  oxidized  ACFC  was  not  related  to  the 
level  of  activation  of  the  precursor  virgin  ACFC. 

As  is  the  case  for  BET  surface  area,  the  micropore  and 
total  pore  volumes  for  the  hydrogen  treated  samples  were 
comparable  but  slightly  lower  than  that  for  the  virgin  sam¬ 
ples,  while  the  acid  treated  samples  depicted  much  lower  mi¬ 
cropore  and  total  pore  volumes  compared  to  the  virgin  sam¬ 
ples.  In  terms  of  ratio  of  micropore  to  total  pore  volume,  the 
acid-treated  samples  depicted  a  reduction  in  microporosity 
for  all  ACFC  activation  levels.  The  hydrogen-treated  sam¬ 
ples  manifested  reduction  in  microporosity,  particularly  for 
ACFC  1  OH  and  ACFC25H. 

For  all  ACFC  grades,  acid  treatment  resulted  in  more 
change  in  mesopore  size  distribution  compared  to  hydro¬ 
gen  treatment.  Acid  treatment  increased  the  mesopore  size. 
Such  change  resulted  in  increasing  the  average  pore  size 


0  0.2  0.4  0.6  o.e  1 

P/Ps 


Fig.  3  Water  vapor  isotherms  at  25  °C  on  Calgon  BPL  and 
ASZM-TEDA  (C2A1) 


Fig.  4  Water  vapor  isotherms  at  25  °C  on  ACnF,  FM10,  and 
FM 10-700 

of  the  acid- treated  ACFC.  On  the  other  hand,  hydrogen 
treatment  doesn’t  seem  to  affect  the  mesopore  size  distri¬ 
bution.  With  respect  to  the  micropore  size  distribution,  the 
impact  of  acid  treatment  was  most  pronounced  for  ACFC20 
and  ACFC25  while  ACFC  10  and  ACFC  15  depicted  smaller 
change  in  both  micropore  and  mesopore  size  distributions. 
For  all  ACFC  grades,  micropore  volume  was  much  larger 
than  the  mesopore  volume  reflecting  the  high  microporosity 
of  ACFC. 
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Table  2  ACFC  composition  as  determined  by  XPS 


Sample  name 

Carbon  (%) 

Oxygen  (%) 

Nitrogen  (%) 

Sulfur  (%) 

ACFC-10V 

94.1 

5.4 

0.3 

0.3 

ACFC- 15  V 

95.3 

3.6 

1.0 

0.1 

ACFC-20V 

94.2 

4.8 

0.6 

0.3 

ACFC-25V 

95.4 

4.3 

0.2 

0.2 

ACFC-10A 

75.1 

23.8 

0.9 

0.2 

ACFC- 15  A 

74.4 

24.3 

1.0 

0.2 

ACFC-20A 

72.6 

25.3 

1.8 

0.3 

ACFC-25A 

76.4 

22.3 

1.2 

Not  available 

ACFC- 1  OH 

98.4 

1.2 

0.4 

Not  available 

ACFC-15H 

98.7 

1.1 

0.2 

Not  available 

ACFC-20H 

97.8 

1.6 

0.6 

Not  available 

ACFC-25H 

97.5 

2.1 

0.4 

Not  available 

Table  3  Characterization  of  ACFC  samples  using  N2  adsorption 


Sample 

BET 

surface  (m2/g) 

Micropore 
volume  (cm3/g) 

Total  pore 
volume  (cm3/g) 

Microporosity  (%) 

Average  micropore 

size  (nm) 

Average  pore 

size  (nm) 

ACFC  10- V 

849 

0.40 

0.40 

99.3 

0.67 

0.67 

ACFC15-V 

1335 

0.62 

0.64 

96.6 

0.76 

0.76 

ACFC20-V 

1566 

0.74 

0.76 

97.5 

0.84 

0.86 

ACFC25-V 

1763 

0.83 

0.84 

98.9 

0.88 

0.88 

ACFC  10-0 

793 

0.33 

0.38 

86.1 

0.70 

0.73 

ACFC  15-0 

994 

0.45 

0.49 

91.9 

0.82 

0.94 

ACFC20-O 

668 

0.29 

0.33 

85.2 

0.82 

1.08 

ACFC25-0 

873 

0.37 

0.44 

85.0 

0.85 

1.11 

ACFC10-H 

926 

0.39 

0.44 

88.4 

0.71 

0.71 

ACFC15-H 

1314 

0.60 

0.62 

95.6 

0.80 

0.80 

ACFC20-H 

1509 

0.70 

0.72 

97.6 

0.86 

0.86 

ACFC25-H 

1739 

0.77 

0.83 

93.1 

0.88 

0.88 

While  the  average  micropore  size  was  less  influenced  by 
the  chemical  treatment,  the  average  pore  size  for  pores  be¬ 
tween  4  and  100  A  was  much  larger  for  acid- treated.  This 
confirms  the  above  suggestion  that  acid  treatment  resulted 
in  micropores  fusing  together  to  form  larger  pores.  On  the 
other  hand,  both  virgin  and  hydrogen-treated  samples  de¬ 
picted  similar  average  micropore  and  total  pore  sizes  indi¬ 
cating  small  influence  of  hydrogen  and  heat  treatment  on 
the  pore  structure  of  ACFC.  The  decrease  in  the  micropore 
volume  and  the  absence  of  prominent  effect  on  the  microp¬ 
ore  size  is  consistent  with  previous  results  reported  by  other 
researchers  (Mangun  et  al.  1999). 

To  summarize,  hydrogen  treatment  has  very  small  im¬ 
pact  on  micropore  and  total  pore  volumes,  average  micro¬ 
pore  size,  BET  surface  area.  On  the  other  hand,  acid  treat¬ 
ment  decreases  the  micropore  and  total  pore  volumes,  mi¬ 
croporosity,  BET  surface  area.  However  acid  treatment  in¬ 


creases  the  average  pore  size  and  doesn’t  have  significant 
effect  on  the  micropore  size.  Hence,  these  results  suggest 
that  the  addition  of  oxygenated  functional  groups  by  treat¬ 
ment  with  HNO3  and  H2SO4  mixture  resulted  in  decreasing 
the  volume  accessible  to  N2.  It  might  be  that  the  functional 
groups  have  blocked  the  pores  of  the  ACFC.  The  large  mi¬ 
cropores  which  are  readily  accessible  to  nitrates  and  sulfate 
ions  might  be  affected  more  by  pore  blockage  than  the  small 
micropores  that  might  be  less  accessible  to  those  ions.  This 
can  justify  the  reduction  in  the  micropore  and  total  pore  vol¬ 
umes.  The  increase  in  the  mesopore  size  might  be  explained 
by  the  action  of  the  concentrated  acid  in  etching  the  car¬ 
bon  structure  and/or  by  fusion  of  large  micropores  and  small 
mesopore  forming  larger  pores. 
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Fig.  5  Water  vapor  isotherms  at  25  °C  on  ACC-5092- 10,  -15,  -20,  and 
-25 


P/Ps 


Fig.  6  Water  vapor  isotherms  at  25  °C  on  acid-treated  ACC-5092- 10, 
-15,  -20,  and  -25 

9.3  Water  adsorption  isotherm  results 

Isotherm  plots  for  the  various  carbons  are  presented  in 
Figs.  3-7. 

The  BPL  isotherm  is  in  good  agreement  with  previously 
published  results  (Huggahalli  and  Fair  1996;  Rudisill  et  al. 
1992).  The  ASZM-TEDA  carbon  is  BPL  with  impregnants 


Fig.  7  Water  vapor  isotherms  at  25  °C  on  hydrogen-treated 
ACC-5092-10,  -15,  -20,  and  -25 

added  to  enhance  adsorption  of  chemical  warfare  agents. 
The  impregnants  add  primary  adsorption  sites  and  reduce 
the  total  pore  volume.  Thus  the  onset  of  water  adsorption 
for  the  ASZM-TEDA  is  earlier,  and  the  total  water  adsorp¬ 
tion  is  reduced  compared  to  the  BPL. 

The  ACnF  exhibits  minimal  hysteresis  and  is  the  most 
hydrophilic  of  all  the  carbons  tested.  Because  the  precursor 
of  this  material  is  PAN,  the  material  likely  has  more  primary 
adsorption  sites  than  the  non-nitrogen  containing  carbons. 
Preliminary  XPS  results  for  ACnF  show  8%  total  nitrogen. 

The  Calgon  Zorflex™  FM10  and  FM 10-700  use  rayon 
as  a  precursor.  Their  curves  most  closely  resemble  those 
of  BPL,  whose  starting  material  also  contains  a  significant 
quantity  of  cellulose  structures. 

The  Kynol™  ACFCs  exhibit  a  strong  type  5  isotherm  ac¬ 
cording  to  IUPAC  classification.  This  is  characterized  by 
low  adsorption  capacity  at  low  RH  followed  by  a  sharp  rise 
and  then  a  plateau  with  little  increase  in  adsorption  capac¬ 
ity  with  increasing  relative  humidity.  The  degree  of  hystere¬ 
sis  increases  with  the  level  of  activation.  The  shape  of  the 
ACFC  isotherms  most  closely  resemble  those  of  pitch-based 
carbons  examined  by  Kaneko  et  al.  (Kimura  et  al.  2004; 
Ohba  and  Kaneko  2007).  This  is  expected,  as  the  pitch  and 
phenolic  resin  carbons  are  chemically  similar,  as  in  the  case 
of  the  BPL  and  Zorflex™.  The  Kynol™  is  a  slightly  more 
homogenous  material,  and  so  the  highly  activated  samples 
exhibit  more  idealized  behavior,  specifically  the  separation 
at  the  bottom  of  the  hysteresis  loop  is  greater. 

The  oxidation  of  the  ACFCs  transformed  the  isotherms 
from  type  5  IUPAC  classification  to  type  4.  The  materials 
were  drastically  modified  in  both  chemical  composition  and 
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Fig.  8  Water  vapor  adsorption 
isotherms  at  25  °C  on  virgin 
ACC-5092- 10,  with  plots  of 
fitted  equations 
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Fig.  9  Water  vapor  adsorption 
isotherms  at  25  °C  on 
acid-treated  ACC-5092- 15,  with 
plots  of  fitted  equations 
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pore  structure  as  noted  in  Table  3,  and  the  isotherms  re¬ 
flect  this.  With  >20%  oxygen,  the  primary  adsorption  sites 
are  plentiful  and  water  adsorption  increases  with  increas¬ 
ing  RH  throughout  the  isotherm.  Although  this  treatment 
yields  increased  water  adsorption  in  high-humidity  situa¬ 
tions,  which  may  be  undesirable,  it  has  utility  in  also  in¬ 


creasing  adsorption  of  other  difficult- to-adsorb  compounds, 
such  as  hydrides  (Mangun  et  al.  1999). 

The  hydrogen-treated  ACFCs  have  increased  hydropho- 
bicity.  The  -10  material  is  most  radically  changed,  going 
from  negligible  hysteresis  to  significant  hysteresis.  As  the 
pore  size  is  relatively  unchanged,  this  is  an  indication  that 
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Fig.  10  Water  vapor  adsorption 
and  desorption  isotherm  at  25  °C 
on  acid-treated  ACC-5092- 15, 
with  plots  of  fitted  equations 


Fig.  11  Water  vapor  adsorption 
isotherms  at  25  °C  on 
hydrogen-treated  ACC-5092-20, 
with  plots  of  fitted  equations 
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the  surface  chemistry  is  key  to  the  transformation.  The  -15 
material  also  is  altered  significantly,  with  the  hysteresis 
loop  being  widened  at  the  base.  The  -20  and  -25  materi¬ 
als  show  a  less  dramatic  general  widening  of  the  hysteresis 
loop. 


When  Figs.  2  and  7  are  compared,  it  is  apparent  that  the 
hydrogen-treated  Kynol™  ACFC  is  more  similar  to  the  ide¬ 
alized  slit-shaped  pore  model  than  previously  studied  car¬ 
bons.  The  hysteresis  loop  is  wider  and  more  squared  at  the 
bottom  of  the  loop.  The  plateau  is  of  a  similar  shape  and 
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Fig.  12  Water  vapor  adsorption 
and  desorption  isotherms  at 
25  °C  on  hydrogen-treated 
ACC-5092-20,  with  plots  of 
fitted  equations 


Fig.  13  Water  vapor  adsorption 
isotherm  at  25  °C  on  Calgon 
BPL,  with  plots  of  fitted 
equations 
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slope.  The  adsorption  loops  commence  adsorption  around 
an  RH  of  60%. 

The  Liu  and  Monson  model  predicts  that  if  there  are 
no  primary  adsorption  sites,  condensation  can  still  oc¬ 


cur  in  the  pore  below  the  saturation  pressure  if  the  pore 
width  is  optimal — about  0.65  nm.  The  ACFCs  have  a 
pore  size  distribution  that  is  not  homogenous,  and  so  there 
are  some  pores  at  this  optimal  size.  The  water  that  con- 
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Fig.  14  Water  vapor  desorption 
isotherm  at  25  °C  on  Calgon 
BPL,  with  plots  of  fitted 
equations 
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Fig.  15  Water  vapor  adsorption 
and  desorption  isotherm  at 
25  °C  on  Calgon  BPL,  with 
plots  of  fitted  equations 


denses  in  the  optimal  pores  then  becomes  a  secondary 
adsorption  site  for  further  water  adsorption.  Thus,  sig¬ 
nificant  water  adsorption  would  generally  commence  at 
a  RH  of  ~60%  in  all  supermicroporous  materials  re¬ 
gardless  of  the  surface  chemistry,  based  on  the  GEMC 
model. 


10  Equilibrium  isotherm  modeling  results 

We  applied  a  least  squares  fitting  algorithm  to  the  equations 
listed  in  Table  1,  and  calculated  the  r2  values  of  the  equa¬ 
tions  against  the  experimental  data.  The  results  are  summa¬ 
rized  in  Table  4. 
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Fig.  16  Water  vapor  adsorption 
and  desorption  isotherm  at 
25  °C  on  ASZM-TEDA,  with 
plots  of  fitted  equations 


Fig.  17  Water  vapor  adsorption 
isotherm  at  25  °C  on  eSpin 
ACnF,  with  plots  of  fitted 
equations 
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Most  equations  seem  to  fit  well,  at  least  numerically,  to 
the  hydrophilic  carbons.  It  can  be  seen  from  Table  4  that 
the  virgin  and  hydrogen-treated  Kynol™  fabrics  have  overall 
worse  correlation  with  these  equations  than  the  other  adsor¬ 
bent  materials.  There  are,  however,  some  equations  that  do 


well  on  both:  the  CMMS-type  equations  and  the  Stoeckli, 
CIMF,  QHR,  Talu-Meunier,  and  Mahle  equations  all  have 
minimum  r 2  values  of  greater  than  99%  for  this  data  set. 
In  fact,  the  largest  r2,  99.99%,  occurs  for  one  of  these  hy¬ 
drophobic  carbons  with  the  full  CMMS  equation. 
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Table  4  The  regression  coefficient  r2  for  each  pair  of  equation  and  experimental  data 
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DR 

DA 

Stoeckli 

DS2 

DS4 

LeVan 

DoDo 

CIMF 

Cossarutto 

ACnF 

95.53% 

99.18% 

99.61% 

98.23% 

99.85% 

99.39% 

99.81% 

99.81% 

99.93% 

ASZM-TEDA  Canister 

98.76% 

98.79% 

99.86% 

99.14% 

99.75% 

99.93% 

99.57% 

99.79% 

99.86% 

BPL 

99.49% 

99.83% 

99.78% 

99.36% 

93.59% 

99.31% 

98.25% 

99.95% 

99.88% 

Zorflex  FM10 

99.65% 

99.70% 

99.90% 

98.86% 

99.37% 

99.42% 

99.83% 

99.88% 

99.88% 

Zorflex  FM 10-700 

99.14% 

99.70% 

99.62% 

99.32% 

99.31% 

99.31% 

99.78% 

99.79% 

99.79% 

ACFC10V 

95.05% 

99.62% 

99.51% 

95.40% 

99.86% 

98.47% 

96.21% 

99.84% 

99.94% 

ACFC15V 

94.75% 

99.60% 

99.49% 

97.91% 

99.46% 

98.41% 

92.69% 

99.79% 

99.90% 

ACFC20V 

94.54% 

99.77% 

99.71% 

98.28% 

99.58% 

98.81% 

91.56% 

99.90% 

99.93% 

ACFC25V 

93.13% 

99.95% 

99.94% 

97.78% 

89.25% 

98.42% 

88.79% 

99.93% 

99.93% 

ACFC10H 

92.82% 

99.89% 

99.97% 

97.47% 

98.79% 

98.02% 

88.37% 

99.94% 

99.94% 

ACFC15H 

94.26% 

99.94% 

99.93% 

98.68% 

93.37% 

98.77% 

88.39% 

99.95% 

99.96% 

ACFC20H 

98.26% 

99.77% 

99.70% 

99.54% 

99.52% 

99.52% 

93.49% 

99.81% 

99.82% 

ACFC25H 

98.15% 

99.73% 

99.69% 

98.10% 

89.93% 

98.13% 

92.09% 

99.80% 

99.81% 

ACFC15A 

99.36% 

99.33% 

99.84% 

99.27% 

98.98% 

99.22% 

99.87% 

99.85% 

99.89% 

Max 

99.65% 

99.95% 

99.97% 

99.54% 

99.86% 

99.93% 

99.87% 

99.95% 

99.96% 

Min 

92.82% 

98.79% 

99.49% 

95.40% 

89.25% 

98.02% 

88.37% 

99.79% 

99.79% 

Lagorsse 

CMMS  full 

XCMMS 

CMMSV  QHR 

Talu 

Mahle 

Max 

Min 

ACnF 

99.12% 

99.73% 

99.85% 

99.75% 

99.58% 

98.61% 

99.87% 

99.93% 

95.53% 

ASZM-TEDA  Canister 

98.71% 

99.45% 

99.87% 

99.49% 

99.46% 

98.92% 

99.76% 

99.93% 

98.71% 

BPL 

98.01% 

99.94% 

99.78% 

99.82% 

99.84% 

99.50% 

99.50% 

99.95% 

93.59% 

Zorflex  FM10 

98.68% 

99.83% 

99.83% 

99.84% 

99.84% 

98.97% 

99.73% 

99.90% 

98.68% 

Zorflex  FM  10-700 

98.77% 

99.77% 

99.68% 

99.72% 

99.76% 

99.06% 

99.49% 

99.79% 

98.77% 

ACFC10V 

98.42% 

99.93% 

99.97% 

99.94% 

99.84% 

98.99% 

99.98% 

99.98% 

95.05% 

ACFC15V 

95.89% 

99.90% 

99.93% 

99.91% 

99.83% 

99.03% 

99.92% 

99.93% 

92.69% 

ACFC20V 

94.57% 

99.95% 

99.95% 

99.96% 

99.93% 

99.36% 

99.93% 

99.96% 

91.56% 

ACFC25V 

92.93% 

99.97% 

99.91% 

99.92% 

99.91% 

99.66% 

99.80% 

99.97% 

88.79% 

ACFC10H 

92.67% 

99.99% 

99.97% 

99.97% 

99.95% 

99.23% 

99.87% 

99.99% 

88.37% 

ACFC15H 

91.55% 

99.89% 

99.88% 

99.90% 

99.95% 

98.98% 

99.72% 

99.96% 

88.39% 

ACFC20H 

93.60% 

99.79% 

99.77% 

99.80% 

99.88% 

98.98% 

99.52% 

99.88% 

93.49% 

ACFC25H 

91.78% 

99.80% 

99.78% 

99.81% 

99.88% 

99.06% 

99.44% 

99.88% 

89.93% 

ACFC15A 

96.94% 

99.70% 

99.90% 

99.69% 

99.36% 

98.98% 

99.75% 

99.90% 

96.94% 

Max 

99.12% 

99.99% 

99.97% 

99.97% 

99.95% 

99.66% 

99.98% 

99.99% 

Min 

91.55% 

99.45% 

99.68% 

99.49% 

99.36% 

98.61% 

99.44% 

88.37% 

Max  and  Min  are  the  maximum  and  minimum  r 2  value  along  each  row  or  column.  The  Max/Max  and  Min/Min  cells  contain  the  maximum  and 
minimum  r2  value  for  the  entire  table 


The  r2  values  can  be  misleading,  though.  It  is  enlighten¬ 
ing  to  look  at  the  graphs  themselves  and  look  for  aberrations. 
The  DR  and  Do  &  Do  equations  can  be  seen  to  have  very 
significant  overshoot  near  RH  =  100%.  This  is  occasionally 
experienced  with  the  CIMF,  and  full  CMMS  as  well.  The 
DR  and  Do  &  Do  equations  also  fail  to  capture  the  steep¬ 
ness  of  the  S-shaped  part  of  an  adsorption  curve.  The  LeVan 
equation  sometimes  returns  a  fit  that  is  not  monotonic,  such 
as  in  ACFC10V.  The  Lagorsse  equation  does  not  fit  very 
well  on  hydrophobic  carbons  such  as  ACFC20H,  and  for 


hydrophilic  carbons  such  as  ACFC15S  often  gives  a  desorp¬ 
tion  branch  that  is  lower  than  the  adsorption  branch,  which 
is  an  aphysical  result.  Several  equations  show  strange  jumps 
sometimes,  such  as  Stoeckli  and  Do  &  Do  and  CIMF  on 
eSpin  near  RH  =  0%,  and  the  DS4  on  BPL. 

It  is  important  to  remember  that  these  correlations  are  for 
these  specific  carbons,  and  in  fact  most  of  the  equations  that 
perform  poorly  here  usually  performed  well  on  the  carbons 
examined  in  the  papers  in  which  these  formulas  were  pro¬ 
posed. 


4^  Springer 


188 


Adsorption  (2007)  13:  173-189 


11  Conclusions 

There  is  a  variety  of  equilibrium  isotherm  models  in  the  lit¬ 
erature.  These  models  were  compared  over  a  diverse  group 
of  activated  carbons.  Although  some  of  these  models  can  fit 
the  isotherms  with  an  r2  greater  than  99%,  none  of  these 
models  are  predictive.  There  are  no  isotherm  models  that 
model  hysteresis  with  any  reasonable  degree  of  accuracy. 

Observing  water  adsorption  isotherms  with  hysteresis 
gives  qualitative  information  about  the  chemistry  and  struc¬ 
ture  of  the  carbons.  Carbons  that  are  similar  chemically 
have  similarly-shaped  isotherms,  e.g.,  carbons  from  cellu- 
losic  materials. 

Novoloid-based  ACFC  with  hydrogen  treatment  is  more 
hydrophobic  and  has  a  more  idealized  isotherm  than  any 
previously  presented  activated  carbon.  This  material  would 
be  ideal  for  high-humidity  adsorption  applications. 

Although  no  real  material  can  be  modeled  as  simply  as  a 
slit-shaped  pore,  the  ACFC  is  an  interesting  material  for  the 
subject  of  further  molecular  adsorption  simulations. 
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